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Abstract 

The Davey-Stewartson equations are used to describe the long time evolution of a three- 
dimensional packets of surface waves. Assuming that the argument functions are quadratic in 
spacial variables, we find in this paper various exact solutions modulo the most known symmetry 
transformations for the Davey-Stewartson equations. 
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1 Introduction 

The Davey and Stewartson [6] (1974) used the method of multiple scales to derive the 
following system of nonlinear partial differential equations 

2iu t + e\U xx + u yy — 2e 2 \u\ 2 u — 2uv = 0, (1.1) 

v xx - e^Vyy + 2(\u\ 2 ) xx ) = (1.2) 

that describe the long time evolution of a three-dimensional packets of surface waves, 
where u is a complex- valued function, v is a real valued function and ei,e2 = ±1. 
The equations are called the Davey-Stewartson I equations if ei = 1, and the Davey- 
Stewartson II equations when E\ = —1. They were used to study the stability of the 
uniform Stokes wave train with respect to small disturbance. The soliton solutions of the 

^000 Mathematical Subject Classification: Primary 35C05, 35Q55; Secondary 37K10. PACS: 4710Fg, 
47.35.Fg. 

2 Research supported by China NSF 10871193 



Davey-Stewartson equations were first studied by Anker and Freeman [2] (1978). Kirby 
and Dalrymple [7] (1983) obtained oblique envelope solutions of the equations in interme- 
diate water depth. Omote [13] (1988) found infinite-dimensional symmetry algebras and 
an infinite number of conserved quantities for the equations. 

Arkadiev, Pogrebkov and Polivanov [3] (1989) studied the solutions of the Davey- 
Stewartson II equations whose singularities form closed lines with string-like behavior. 
They [4] (1989) also applied the inverse sacttering transform method to the Davey- 
Stewartson II equations. We refer [1] for more information on the Davey-Stewartson 
equations in terms of inverse scattering and integrable systems. Gilson and Nimmo [8] 
(1991) found dromion solutions and Malanyuk [11, 12] (1991, 1994) obtained finite-gap 
solutions of the equations, van de Linden [15] (1992) studied the solutions under a certain 
boundary condition. Clarkson and Hood [5] (1994) obtained certain symmetry reductions 
of the equations to ordinary differential equations with no intervening steps and provided 
new exact solutions which are not obtainable by the Lie group approach. Guil and Manas 
[9] (1995) found certains solutions of the Davey-Stewartson I equations by deforming 
dromion. Manas and Santini [10] (1997) studied a large class of solutions of the Davey- 
Stewartson II equations by a Wronskian scheme. It is obvious that the some of above 
solutions are equivalent to each other under the known symmetric transformations. It is 
time to study solutions of the Davey-Stewartson equations modulo the known symmetric 
transformations. 

We observed in [17] that the argument functions of all the solutions of the two- 
dimensional cubic nonlinear Schrodinger equation in [15] are quadratic in spacial variables. 
This motivated us [17] to introduce a quadratic- argument approach to study exact solu- 
tions of the two-dimensional cubic nonlinear Schrodinger equation and the coupled two- 
dimensional cubic nonlinear Schrodinger equations modulo the known symmetry trans- 
formations. Indeed, our solution sets are most complete among the ones whose argument 
functions are quadratic in spacial variables. In this paper, we use the quadratic-argument 
approach to study exact solutions of the Davey-Stewartson equations modulo the most 
known symmetry transformations. 

For convenience, we always assume that all the involved partial derivatives of related 
functions always exist and we can change orders of taking partial derivatives. We also 
use prime ' to denote the derivative of any one-variable function. The known symmetry 
transformations we are concerned with are: 



Ti(m) = e -W x+fy **H{t, x + a,y + P), 
T\(y) = v(t, x + a,y + (3) + e\ol' x + f3''y — — ^—^ 




(1.3) 



(1.4) 



T 2 (u) = 6 -1 u(6 ? t, bx, by), T 2 (v) = b~ 2 v(b 2 t, bx, by); (1.5) 

where a, (3, 7 are functions oft and 6 is a nonzero real constant. The above transformations 
transform one solution of (1.1) and (1.2) into another solution. In particular, applying 
the transformation Ti to any known solution would yield solutions with three additional 
parameter functions. 

2 Exact Solutions 

Write 

u = Z(t,x,y)e i *<t>*>v), (2.1) 
where £ and <fi are real functions in t,x,y. Note 

u t = (& + i^t)^, u x = (£, + i^e* = (£„ + i^e* (2.2) 

«,,■ = (U - O-r + iV>l, o, + ^ xx ))e i4 ', u yy = (i yy - + *(2£A + £0 TO ))e^. (2.3) 
Then (1.1) is equivalent to 

2z£ t - + - icfi + i(2^0 x + & xx )) 

Hyy - + i(^ y (f>y + & m ) ~ 2e ^ 3 ~ 2 ^ = °> ( 2 - 4 ) 

equivalently, 

26 + 2(e 1 £ E a; + i v 4> y ) + f (e^ + W ) = 0, (2.5) 

£(20 t + t l( f>l + 2 ) - e^xx - Zyy + 2e 2 £ 3 + 2£v = 0. (2.6) 
Moreover, (1.2) becomes 

v xx -e 1 (v yy + 2(f) xx ) =0. (2.7) 

In order to solve the above system of partial differential equations, we assume in this 
section that 

4> = a 'x 2 + (3'y 2 (2.8) 
for some functions a and (3 of t. Then (2.5) becomes 

& + 2(e 1 a , z6 E + /?'y4) + ( £l a' + (3% = 0. (2.9) 

Thus 

£ = e - eia_ ^(n7i, C7 2 ), G7i = e" 2eiQ x, ro 2 = e~ 2/3 2/, (2.10) 
where i? is some two-variable function in w\ and w 2 . Moreover, (2.6) is equivalent to 

2((o" + 2e 1 (o / ) V + ([3" + 2{(3>) 2 )y 2 ¥ ~ 

- e - 4/3 ^ ro2ro2 + 2e 2 e~ 2eia - 2 ^ 3 + 2-&v = 0. (2.11) 



Case 1. d — aw\ + bzu 2 + c for a,b,c G R. 

In this case, The above equation is equivalent to 

(a" + 2 ei (a') V + + 2(/3') V + € 2 e- 2 ^ 2 f (ae~ 2eia x + fre" 2 ^ + c) 2 + v = 0. (2.12) 
So 

v = -[«" + 2e 1 (a') 2 + e 2 aV 6eia ~ 2 V - [/?" + 2{(3'f + e 2 6 2 e - 2eiQ - 6 V 

-2a6e 2 e- 4eia - 4/3 x ?/ - t 2 ce- 2eia - w {2ae~ 2 ^ a x + 2be~ 2li y + c). (2.13) 

Moreover, (2.7) is equivalent to 

a" + 2 ei («') 2 + e 2 a 2 e"^ a ^ = e 1 [p" + 2((3') 2 + e 2 b 2 e~ 2 ^ a -^\. (2.14) 

To solve the above equation, we write 

a = ^s + e l (3. (2.15) 

Then the above equation becomes 

3" + 2e 1 (S / ) 2 + 43'/?' + e 2 (aV 6eiS - ei feV 2ei3 ) e - 8/3 = 0. (2.16) 
Assuming ab ^ and ei = 1, we have the following simple solution 

3f = i(ln|a| — In j&|), (2.17) 

where ft is any function of t. Suppose a = b = 0. Then we take 

cy' , cy i , c> 

/?' = -—- ^ =► = -- In 3' - ^ (2.18) 

^ 43/ 2 M 4 2 V ; 

modulo the transformation in (1.5), where 3 is arbitrary increasing function of t. Then 

a = (3')~ £l/4 e 3/2 , p = (s /)-/4 e - ei a/2. ( 2 .19) 

By (2.8), (2.10) and (2.13), we have: 

Theorem 2.1. Let S &e any increasing function of t and let a,b,c G R with ab ^ 
0. Suppose that (3 is any function of t. We have the following solutions of the Davey- 
Stewartson equations (1.1) and (1.2): 

u = cexp(- ei (^)" ei/4 e 3/2 -(^)~ /4 e^ l3/2 ) 

(2(3') 2 - ei 3'V - (2ei(3') 2 + 3'V . 

x exp ^ ^ U ^ J —i, (2.20) 



ei(S'''- (S'Q 2 ) -2(S') 2 S'' 2 S"'- (^) 2 + 2e 1 (^ / ) 2 S // 2 
V ~ 4(3') 2 ^ + 4(3') 2 ^ 

-e 2 c 2 exp(-2 ei (S')^ l/4 e 3/2 - 2(^')~/ 4 e-^ 2 ); (2.21) 

u = ^'(x 2 + y 2 )i-(2/3 + (ln|a|-ln| 6 |)/2) (ae ln| 6 |-ln|a|-2^ + ^-2/^ + ^ (2>22) 

-e 2 ce- 4/3 (2a- 1 6 2 e- 2 ' 3 x + 2|a|- 1 6|b| e - 2/3 ?/ + c). (2.23) 



Remark 2.2. Applying the symmetry transformation T\ in (1.3) and (1.4) to the 
above solutions, we get the following solutions with additional three parameter functions 
ai,/?i,7i of t: 



/ 7 7V^ 1 "17 7 i ^ J\y i M17 / nl \ ■ 

ii : cexp( — eia x x - fay - 71 ] % 



(2(3') 2 - e^")(x + a 1 ) 2 _^_(2e 1 ^') 2 + %")(y + /30 2 

43' 

x exp(- ei (^)^ l/4 e 3/2 - (S')~ /4 e~ ei3/2 ), , (2.24) 



e^Q"' -(&')*) -2(&)*Q", a ^ - (^) 2 + 2 £l (^) 2 ^ 2 „ 



+ fa' y _ e2C 2 exp(-2e 1 (3 / )- £l/4 e 3 / 2 - 2^'y/ 4 e -^ 2 ) - !l^il!±Mj! + ^ ; 



(2.25) 



if ei = 1, 



it = e 



[^'(x 2 +j/ 2 )-eia' 1 x+^j/+7i]i-(2^+(ln |a|-ln |6|)/2) 



x (ae ln l b l- ln l a '- 2/3 (x + a ± ) + be- 2 \y + A) + c), (2.26) 



-[fa' + 2(/3') 2 + e 2 \a\^\b\ 3 e-^]((x + Ql ) 2 + (y + fa) 2 ) + e^x + 
-e 2 ce" 4/3 (2a- 1 6 2 e- 2 ' 3 (x + «i) + 2|a|- 1 6|fe| e - 2/3 ( 2 , + ft) + c) 

-2a- 1 6 3 6 2 e- 8 ^(x + ai )(i/ + A) - fl^tMl! + ^ {2 .27) 



Case 2. t\a = fa 
We set 

Ci(s) = sinhs, r]i(s) = coshs, C-i( s ) — sins, r)-i(s) = coss. (2.28) 



Let £ and £\ be two fixed real constants. We denote 

m = ( £1 (£)m 1 + r 1ei (£)m 2 + t 1 . (2.29) 

Moreover, we assume 

= v{w). (2.26) 

Observe 

(d 2 x -e 1 d')(v(zu)) = -e 1 v"(zu) (2.30) 

and 

{e l dl + dl){v{ui)) = {r ltl {2£)y{w). (2.31) 
By (2.7), we further assume 

v = 7 ( ei x 2 + y 2 ) + e-^(c + £ ( V) (2.32) 

for a function 7 of t and a real constant c. So (2.7) naturally holds. 
Now (2.11) becomes 

2(/3" + 2(/3') 2 + 7 )M 2 + y> + e~ 4 V - fciW + 2(e 2 + Cj W)^ 3 ) = 0. (2.33) 
For simplicity, we take 

7 = - 2(/3') 2 . (2.34) 

So (2.29) is equivalent to 

cv - VeA^V + 2(e 2 + CiW)^ = 0. (2.35) 
First we have simple solution 



Recall 

(tan s)'' = 2(tan 3 s + tans), (sec s)'' = 2 sec 3 s — sec s, (2.37) 

(coths)" = 2(coth 3 s - coths), (csch s)" = 2csch 3 s + csch s. (2.38) 
Denote Jacobi elliptic functions 

sns = sn(s|m), cns = cn(s|m), dns = dn(s|m), (2.39) 

where m is the elliptic modulus (e.g., cf. [16]). Then 

(sn s)" = 2m 2 sn 3 s - (1 + m 2 )sn s, (2.40) 



(cn s)" = -2m 2 cn 2 s + (2m 2 - l)cn s, (2.41) 

(dns)" = -2dn 3 s + (2-m 2 )dns. (2.42) 

Moreover, 

lim sn s = tanh s, lim cn s = lim dn s = sech s. (2.43) 

m—*l m—*l m—*l 

Comparing (2.35) with the equations in (2.37), (2.38) and (2.39)-(2.42), we have the 
following solutions: 



u= \l € ?+ ( ?\e ) secro ' o=-r kl (2£); (2.45) 



y= ^ e 2 + CMf) C ° thCT ' c =- 2 ^( 2 ^); ( 2 - 46 ) 



* = m \I ^^Ji) sn w ' c = -(1 + m 2 )r ?£1 (2£); (2.48) 



,/= '"*/ l £l ^ Cnro ' c=(2m 2 -l)r ?ei (2£), (2.49) 



^ l(2 ^ dnro, c= (2-m 2 )r/ ei (2f). (2.50) 



^2+ cm 

By (2.8), (2.10) and (2.32), we have: 

Theorem 2.3. Let £,£i,m e R. Suppose that (3 is any function oft. Take the nota- 
tions in (2.28). The fallowings are solutions of the Davey-Stewartson equations (1.1) and 
(1.2) (where the solution exists only when the expression makes sense as real function) : 



e 



»fci(2*) 



« = kHt^TX e- 2/3+(ei ^ )/3 ' l tan(e- 2 ^(C 1 + r? ei (%) + 4), (2.53) 



e 2 + C, 



-(/?" + 2(/?) 2 )( ei x 2 + y 2 ); 



-09" + 2(/3') 2 )(e 1 x 2 + y 2 ); 



« = J^TT^k e-^^)^ coth(e-^(C 1 + ^)y) + 4), 
e 2 + C 2 (^) 

- 1 ' ( QW^M coth 2 (e- 2 ^(C ei + ifc^v) + €i) - 2^(2*) 



e 2 + C 2 W 
-(/?" + 2(/?') 2 )(^ 2 + y 2 ); 



e " 4 " ( t + c 2 l( ffl csch V^kO* + ^ (%) + 4) + 1^(2*)) 
-(/?" + 2(/3') 2 )(e 1 x 2 + y 2 ); 



e~ AP ((2m 2 - l)r/ ei (2£) - ^^pcn V^C. + r, ei (£)y) + *)) 
-(/?" + 2(/3') 2 )( e ix 2 + y 2 ); 



e-^^^dn (e^(( ei (£)x + ^ (£)y) + 4), 
e2 + C 2 W 



e-^ ( (2 - m 2 )^(2£) - ^^^ dn 2 (e^(( ei (£)x + r/ ei (£)y) + t x ) 



t2+ cm 

-(/?" + 2(/3') 2 )(e 1 x 2 + y 2 ); 



Remark 2.4. Applying the symmetry transformation T\ in (1.3) and (1.4) to the 
above solutions, we can obtain solutions with additional three parameter functions. For 
instance, we get the following solutions with additional three parameter functions cui, fli, 71 
of t from the above first two solutions: 



e -2^+l(e 1 (x+a i y 2 + (y+/3 1 ) 2 )^~e ia ' 1 x-/3[y-y l }i / ^ ^f) 

u = e-^(ue)(x + ai ) + Vei (e)(y + A)) + 4 V e 2 + Cpj ' (2 ' 67) 

t; = -(/?" + 2(/30 2 )(6 1 (x + ai ) 2 + (y + A) 2 ) + ciafo + fly - ei(a ' l)2 + ^ + 71 



& + Q(t))(e- 2f3 ((eM(x + «i) + TkMiv + A)) + 4) 5 



+ ~ I 7^777777^777 77/77^4^ 7 77777 I ^ I ^ (2-68) 



■u = 



»fci(2*) 



2 / 9+[(ei(a;+ai) 2 +(y+/3i) 2 )/3'-eia' ia :-/3( J /-7i]j 



x tan(e- 2/3 (C 1 W(x + a ± ) + ^(^(y + A)) + (2.69) 

« = e" 4/3 ( C ^2 l( ( ^ ) tan 2 (e- 2 ^(C 1 (£)(x + a x ) + V M(v + A)) + *i) + 2^(2*)) + ^ 
-09" + 2(/3') 2 )(ei0r + ai ) 2 + (v + A) 2 ) + ^['x + (3['y - !iMj!±Mj! . ( 2 .70) 
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